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Abstract 



We study the algebra Sp(n,R) of the symplectic model, in particular for the 
cases n=l,2,3, in a new way. Starting from the Poisson-bracket realization 
we derive a set of partial differential equations for the generators as functions 
of classical canonical variables. We obtain a solution to these equations that 
represents the classical limit of a boson mapping of the algebra. We show 
further that this mapping plays a fundamental role in the collective descrip- 
tion of many-fermion systems whose Hamiltonian may be approximated by 
polynomials in the associated algebra, as is done in the simplest versions of 
the symplectic model. The relationship to the collective dynamics is formu- 
lated as a theorem that associates the mapping with an exact solution of 
the time-dependent Hartree approximation. This solution determines a de- 
coupled classical symplectic manifold, thus satisfying the criteria that define 
an exactly solvable model in the theory of large amplitude collective motion. 
The models thus obtained also provide a test of methods for constructing an 
approximately decoupled manifold in fully realistic cases. We show that an 
algorithm developed in one of our earlier works reproduces the main results 
of the theorem. 
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I. INTRODUCTION 



We have been engaged for a decade in an effort to formulate a theory of large amplitude 
collective motion with the special aim of applying it to nuclear physics. The theory has both 
a classical and a quantum dimension. The classical aspect has been most fully developed 
and described in a review [Q. The quantum aspect is presently in a stage of vigorous 
development @J^, supplanting the early work on this part of the theory [H,^. At the same 
time a program of applications to problems of nuclear structure has been undertaken l^j-Q . 

Early in the latter work, we became aware of a paucity of solvable models with some 
physical content. The usefulness of such models is that they provide a testing ground for the 
algorithms that would later be applied to more realistic models. For our initial investigation 
we selected a well-known model of monopole vibrations |]TD[, exactly solvable because the 
Hamiltonian is a polynomial in the generators of the algebra Sp(l,R) (or SU(1,1)). We stud- 
ied this model in two ways. First, by utilizing the classical limit of the algebra, we were able 
to produce an exact solution of the time-dependent Hartree equations (derived previously by 
rather less transparent techniques |]ll|) and by means of this solution a decoupled collective 
Hamiltonian for the monopole vibration. Second, and more important, we could check if the 
same Hamiltonian emerged from the application of the theory of large amplitude collective 
motion. The monopole model provided us with an apparently ideal test of the soundness of 
our algorithms. This test failed initially, forcing us eventually to recognize and correct an 
incompleteness in our previous theory. 

The first goal of this paper is to show that the method developed for the algebra of 
Sp(l,R) can be extended to the algebra of Sp(n,R). In particular, we work out fully the 
cases n=2 and n=3. The former leads to a Hamiltonian with three collective coordinates, 
describing the interaction of a monopole degree of freedom with a quadrupole tensor in 
two spatial coordinates and is thus only of interest as a toy model. On the other hand 
n=3 leads to a Hamiltonian with six degrees of freedom describing both a monopole and a 
three-dimensional quadrupole. This defines the model as not only one of physical interest 



per se, but also because of its connection with the symplectic |jT2| and pseudo-symplectic 
IB| models. The latter, in particular, provides a possibly useful truncation scheme for 
shell-model calculations for other than the lightest deformed nuclei. 

Though not one of the aims of the present paper, this identification will allow us on 
a future occasion to compare the quantum consequences of the collective Hamiltonian to 
be derived in this paper with the results of an exact diagonalization carried out for the 
original many body Hamiltonian. In fact such a comparison has broader implications than 
the accuracy achieved for the special Hamiltonian considered, since it has been demonstated 
that Hamiltonians consisting of suitably chosen polynomials in the generators can give a 
rather precise fit to the low-energy spectra and other properties of even deformed nuclei 



The second aim of this paper is to demonstrate that the extended algorithm formulated 
in connection with the monopole model also provides correct results for the generalized 
models. 

The presentation is organized as follows: In Sec. 2 we give a brief summary of the prop- 
erties of the algebra of Sp(n,R) needed in the ensueing development as well as a discussion 
of the model Hamiltonians to be studied. In the following three sections we then study 
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separately the cases n=l,2,3. We describe in each instance the mapping of the algebra onto 
a classical symplectic manifold, which is tied to the existence of a manifold of solutions of 
the time-dependent Hartree equation and an associated decoupled collective Hamiltonian. 
We then show how the same collective Hamiltonian can be derived from our theory of large- 
amplitude collective motion. The material for the monopole case is a rearrangement with 
different emphasis of results presented previously All other results are new. In Sec. 6, 
we make suggestions for further work, involving both applications and extensions of the re- 
sults of this paper, as well as the study of possible connections with previous research. Two 
appendices, A and C contain important details of the calculations that would impede the 
flow of the argument in the main text. In Appendix B we review briefly the generalization 
of the theory suggested by our previous study of the monopole model and indicate how it 
applies to the more general cases. 



The defining algebra of the group Sp(n,R) is given most simply in terms of n boson 
pairs, Oj, a*, z = 1- ■ - n, where the subscript identifies the destruction and the superscript 
the creation operators. One standard set of generators is composed of three distinct bilinear 
forms in these operators. 



An equivalent set of generators, more useful for the purposes of this paper is given in 
terms of single-particle coordinate and momentum operators, Xi and pj, respectively, related 
to the boson operators in the standard way, namely. 



II. ALGEBRAIC AND DYNAMICAL PRELIMINARIES 



A. The algebra of Sp(n,R) 



B] 




(2.1) 
(2.2) 
(2.3) 




(2.4) 
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(2.5) 



The alternative set of generators takes the form 



Q ij ^i^j} 
■^ij PiPji 




(2.6) 
(2.7) 
(2.8) 
(2.9) 



satisfying the set of commutation relations. 
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[Qij, Kki] 


= i{SjkSii + SjiSik + SikSji + SiiSjk) 






+i{6jkLii + 6jiLik + SikLji + 6iiLjk), 


(2.10) 


[Qij, Lki] 


= \i{.5ljQik — SjkQil + ^ilQjk — ^ikQjl), 


(2.11) 


[Qij, Ski] 


= \i{5ijQik + SjkQii + SiiQjk + SikQji), 


(2.12) 


[Kij, Lki] 


= \i{5iiKjk + 5jiKik - SikKji - 6jkKii), 


(2.13) 


[Kij, Ski] 


= -\i{5iiKjk + 5jiKik + 5ikKji + SjkKn), 


(2.14) 


[Lij, Ski] 


= \i{5ikSji — SjiSik + SiiSjk — SjkSii), 


(2.15) 


[Lij, Lki] 


= \i{^ikLji + 5jiLik — SiiLjk — 6jkLii), 


(2.16) 


[Sij, Ski] 


= —\i{.5ikLji + SjiLik + 5iiLjk + 5jkLii). 


(2.17) 



In terms of these generators, the second order Casimir invariant is given by the expression 

C2 = \{Qiji L-ij} + SijSij LijLij. (2.18) 



B. Dynamical preliminaries 

In the body of the paper we shall find it convenient to assign a different meaning to the 
upper-case symbols used above to denote the various generators. We shall therefore refer to 
these quantities in the text below by the corresponding lower-case letters, Qij — >• qij, etc. We 
shall be concerned with a many-body realization of the algebra, for which we shall use the 
second-quantized formalism. Thus, if ip{yi), '?/'^(x') satisfy the anticommutation relations, 

Mx),^t(x')} = 5(x_x')^ (2.19) 

where x stands for the n-dimensional vector (xi ■ ■ ■ Xn) and any additional intrinsic variables 
such as spin and isospin that are included in the dynamics, then the operators Qij, Kij, etc., 
where, for example, 

Qij = J d^ij^{^)qijij{^), (2.20) 

are the set of generators of Sp(n,R) relevant to the fermion many-body problem. Below we 
shall also suppress the boldface for the quantity x. 

We shall be concerned further with studying the classical limit of Hamiltonian operators 
belonging to the enveloping algebra of Sp(n,R) for n=l,2,3, i.e., Hamiltonians that are 
polynomials in the generators of these algebras. Consider the one-dimensional case, where 
we have just three generators, Qu, Ku, and Su, that we rename by dropping the subscripts. 
We then study a Hamiltonian, 

H = lik + Q) + lKQQ, (2.21) 

that is the sum of an harmonic oscillator part and of a "monopole-monopole" interaction. 
The considerations to be applied to this Hamiltonian can be generalized to one in which 
we add any polynomial in the operator Q, but we shall quote detailed results only for the 
Hamiltonian ( |2.21|) . 
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Consider an arbitrary Slater determinant describing fermions. We shall be interested 
in the expectation value of (|2.21 ) in this state, evaluated in Hartree approximation, i. e. 



to the leading order in A^, the number of particles. Using angular brackets to denote this 
average, we have in this approximation 

(H) = Hc = \{K + Q) + \kQ\ (2.22) 

where, e. g. , 

g = (g), (2.23) 

i. e. , the Hartree average of any generator will be denoted by the same symbol without hat. 
As a further instance, the Hartree approximation to the Casimir invariant of Sp(l,R) is 

(C2) = C2 = -QK + 5l (2.24) 

The most important idea that informs the next section is that we consider Q to be a 
classical collective coordinate and associate it with a corresponding canonical momentum, 
P. The Hartree approximation to the generators then defines a classical limit of the algebra 
in which the commutators are replaced by Poisson brackets. These are a set of partial 
differential equations for the classical generators, Q^K^ and S*, that determine the latter 
as functions of Q and P. This mapping will not only determine the Hartree average as 
a function of Q and P, but will allow us to associate this classical boson mapping with 
a solution of the time-dependent Hartree equation. The concepts mentioned here will be 
rendered precise in the next section and then extended to two and three dimensions in the 
following sections. 



III. THE MONOPOLE MODEL ASSOCIATED WITH THE ALGEBRA OF 

SP(1,R) 

A. Classical mapping and Hartree solution 

As stated in the introduction, the results to be presented in this section have appeared 
in our previous work P|. What we aim for here is a more systematic presentation with 
enhanced emphasis on the significance of the results. We do this by collecting the results 
into a theorem. 

Theorem: For a class of many-Fermion Hamiltonians belonging to the enveloping algebra of 
Sp(l,R), of which Eq. (|2.21|) is a prototype, there exists a two parameter family of Slater 
determinants, defined by density matrices, p[xx'\Q^ P), that describe states belonging to an 
irreducible representation. These states 

(i) induce a mapping of Sp(l,R) onto a symplectic manifold (Q, P) in which the generators, 
(5, are mapped as classical dynamical variables, i.e., 

G'^G(Q,P), (3.1) 

and in particular. 
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Q^Q. (3.2) 

The Hartree expectation value of any member of the enveloping algebra is thereby also 
mapped. 

(ii) The associated density matrices have the diagonal form 

pix,x'\Q,P) = Y.M^\Q,P)rhix'\Q,P), (3.3) 

h 

where the sum is over the occupied orbitals, h. Furthermore tphi^lQ, P) can be written 
as a product 

M^\Q, P) = exp{iPx^)M^, Q). (3.4) 
and (ph is the solution of a constrained Hartree equation 

eh(f>h = (^ - Aa;^)0/i, (3.5) 

n = ^{p^ + x^) + kQx"^, (3.6) 

_ dHc{Q,P = 0) _ dV 

^ - dQ = dQ- ^^-^^ 

Here Ti. is the Hartree Hamiltonian and He is the classical collective Hamiltonian defined in 
( p.22| ). These results identify the density matrix ( |3.3| ) as a solution of the time-dependent 
Hartree equation 

ip = [n,p]. (3.8) 

We turn to the proof of part (i) of the theorem. The Hartree average of a product of two 
generators, Gi and G2, 

{G,G2) = G,{Q,P)G2{Q^P) (3.9) 

can be identified as the leading term in the convolution of these two operators under a 
Wigner transform with respect to the collective variables Q, P. For a consistent evaluation 
of commutators, however, we need the next term in the convolution, 

{[Gi,G2]) ^i[Gi{Q,P),G2{Q.P)]pB, (3.10) 

where the Poisson bracket is to be evaluated with respect to the single canonical pair (Q,P). 
From Eqs. (|2.10|) - (|2.17|) , we thereby obtain the Poisson bracket algebra 

[g,5] = — = 2g, (3.11) 

[Q.K] = ^ = AS, (3.12) 
[K,S] = -2K. (3.13) 

Treated in turn, these differential equations yield the relations 
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S = 2QP, (3.14) 
K = 4QP^ + X{Q), with (3.15) 

0| = -xiQy (3.16) 

Thus 

xiQ) = C/Q. (3.17) 

The constant C can be obtained by evaluating the second-order Casimir invariant in the 
Hartree approximation. We find 

(Cs) = -C = -ArV4, (3.18) 

where the pair of equahties express the results of two separate procedures. On the one hand 
the value —C is obtained by direct substitution of the mapped generators. On the other 
hand the specific value is obtained by calculating the Casimir invariant for a simple 

state in the irreducible representation, as explained in Appendix A. Finally, we record that 
the collective Hamiltonian, which is the Hartree value of the many-particle Hamiltonian, 
maps to 

He = 2QP^ + + IkQ^ + N^SQ, 

= 2QP^ + V{Q), (3.19) 

and we have also displayed the collective potential energy. Note, however, that the singular 
term in the latter, that has a tantalizing resemblance to the scalar Berry potential [Q, 
originates in the many-particle kinetic energy. 

We turn next to part (ii) of the theorem, which provides a construction of the manifold 
of density matrices associated with the classical mapping just given. We show first that the 
single-particle wave functions of which the density matrix is composed are of the form 

M^\Q, P) = exp{iS{x, P))Mx, Q), (3.20) 

i. e, the dependence on the collective momentum is independent of the orbit label. To see 
this, let us calculate the time derivative of the classical variable Q in two ways, directly from 
the classical equations of motion, 

Q = ^ = ^PQ^ (3-21) 

and by actually evaluating the Hartree approximation of the quantum equations of motion, 

Q = -i{[Q,H]) = 2{S) 
f dS 

= 2 / dxx—p{x,x). (3.22) 
J ox 



In obtaining this last result, we have assumed that the orbitals (ph are real. Equations (|3.21 
and (|3.22| ) yield the solution 
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S = Px^. (3.23) 
To find the orbitals (ph introduced in Eq. (|3.20|) , we back up a notch by defining 

= exp{-t f dt' eh[Q{t')]}iJh. (3.24) 
Jo 

The extra adiabatic phase in ( |3.24| ) drops out of the density matrix, but must be included 
in order that ip satisfy the time-dependent Hartree equation, 

= Tiij. (3.25) 

Taking into account both the exphcit and the implicit time dependence contained in ( p. 251 ), 
the latter is replaced by the equation 

ehiph + iP^ + iQ^ = T-t^h- (3.26) 

Substituting the classical equations of motion for the time derivatives and inserting the form 
of iphi the resulting equation has terms of zero, first, and second order in the momentum 
P. The terms of second order are found to cancel, whereas the terms of zero and first order 
give, respectively, the equations 

eh0h = {n- \x^)(t)h, (3.27) 
4Q^ = -0-2.-. (3.28) 

We deal first with ( p.27|) . By means of the defintions 

7Y=i(p2 + ^V), (3.29) 

6^2 = l + 2(Kg-A), (3.30) 

( p. 27] ) becomes 

eh(t)h = 'ficph, (3.31) 

with the normalized solution, in terms of solutions 01^'^°'* for the simple harmonic oscillator 
with unit mass and unit frequency, 

0,(x)=cui/Vf°^(v^x). (3.32) 

It remains only to verify (|3.28|) . This can be done, using the explicit form of u^, derived 
from Eq. ( p.3(J| ) by substituting the value of A = (dV/dQ) from Eq. ( |3.19| ), namely, 

uj^ = C/Q\ (3.33) 

This completes the proof of the theorem. 
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B. Application of the theory of large amplitude collective motion 



It is not the purpose of this section to review yet again the theory of large amphtude 
collective motion. It has been described exhaustively in the works referred to in the in- 
troduction. In so far as the present discussion is concerned the principal aim is to find 
the collective Hamiltonian from this theory. The first step of the procedure is to assume 
that the collective variable is the expectation value of a one-body operator. For the model 
Hamiltonian, Eq. ( p.21| ), with a separable interaction, the logical starting choice is always 



determined by the ingredients of that interaction. In the present instance, the choice is 
the Hartree expectation value of Q, or in other words it is determined by the mapping 
Q Q. This leads automatically to the constrained Hartree equation solved in part (ii) of 
the theorem. 

The next step is to compute the many-particle Hartree energy associated with the filling 
of these orbits with particles. We outline the calculation, in which this Hartree energy is 
identified as the potential energy, V{Q), of the system. We have first 

V{Q) = {H) = tr[i(p2 ^ ^2)^] ^ i^g2_ ^3 34) 

This equation can be transformed into a more useful version, as follows: We first calculate 
the sum of the single-particle energies, 

^e, = tr(7^p) 

h 

= tr[i(/ + x^)p] + kQ^ - Q{dV/dQ). (3.35) 
We then invoke the virial theorem, in the form 

J2eh = uhi{x^p) =u^Q, (3.36) 

h 

which allows us to replace ( pj.35D by the expression 

dQ 



tT[^ip^ + x^)p]=iu^Q-KQ^ + Q—. (3.37) 



When the latter is now substituted into ( |3.34| ), the result is a differential equation for the 
potential energy, 

ViQ) = Q + ^KQ'-Q^, (3.38) 

of which the solution is 

V{Q) = IQ + ^f^Q^ + C/Q. (3.39) 



The last term is the solution of the homogeneous version of (|3.38 ), depending on the same 
constant that appeared in Eq. (|3.171 ). 



The next stage in this procedure is to calculate the collective mass, B, defined by writing 
the collective kinetic energy in the form \BP^. We shall calculate B in two (equivalent) 
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ways. The first, that we have not used in our previous work is connected to the work 
done in proving the theorem of the previous section. Its deceptive simphcity carries special 
instruction to which we shall return below. The procedure is to make use once more of 
the two values of Q found in Eqs. ( p.21|) and (|3.22|) , without assuming that we know the 
coefficient of P^. This yields the equality 

BP = 2{S) = 25. (3.40) 

Since we have already solved for the value of 5* from the PB form of the algebra, we may 
substitute from (|3.14|) into ( p.40|) and verify the value B = AQ found previously. (We are 



not depending on the theorem. If we didn't know 5 as a classical dynamical variable, we 
would calculate it now.) . According to the general theory, the fact that B depends only on 
Q verifies that we are dealing with an exactly decoupled mode. This means that we have 
chosen the collective coordinate correctly, i.e., self consistently. 

We shall next attempt to convince the reader that we have not engaged in sleight of 
hand. We can do this by complicating the model in what appears to be a minimal way, but 
taking it outside of the algebra Sp(l,R). We do this by adding to H a. term ^k^QI, where 
Q4 is the many body version of the operater = x^. Even if we crank only with x^, we 
now find that the potential energy is a function both of Q and of Q4, so that the minimum 
number of collective variables that we can introduce is two. We must therefore apply the 
procedure followed above to calculate the mass to the two-dimensional domain with kinetic 
energy, T, given by 

T = ^B^^^P^ + B^^PPi + \B^^PiPi. (3.41) 

We shall not carry this through in detail, but, in fact, it is a perfectly straightforward exercise 
and yields the mass matrix 

5°° = 4g, (3.42) 
5°^ = 8Q4, (3.43) 
B^^ = 16Qe, (3.44) 

where Qq is the Hartree average of x^. The point is that in the classical collective Hamil- 
tonian, Qq is independent of the two variables introduced previously. Thus the effort to 
decouple a finite symplectic manifold has failed. There is no cure for this ailment. Thus 
we might try to introduce Qq as an additional collective variable by cranking on and 
thus bringing it into the potential energy. By extension with what has been found above, 
however, the augmented kinetic energy will bring in still higher powers of x. 

The method of calculating the collective mass developed in our previous work that 
gives the same result as above is described in Appendix B. 



IV. MONOPOLE PLUS TWO-DIMENSIONAL MODEL ASSOCIATED WITH 

THE ALGEBRA OF SP(2,R) 
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A. Classical mapping and Hartree solution 



We turn to the problem of stating and proving a theorem that generahzes the results 
of the previous section. We first specify the model to be studied. Here it is convenient to 
replace the coordinate generators, Qij by a monopole Qq, and a two-dimensional quadrupole, 
Qi, i = 1,2, where 



Qo = y + (4.1) 

Qi = Jij^ix^ixl-xD^Pix), (4.2) 
Q2= [ ^\x)xiX2^{x). (4.3) 



The Hartree expectation values of these operators will be interpreted as their classical maps. 
We consider one further definition, namely, 

ko = U^n + K22). (4.4) 

As model Hamiltonian we choose 

H = liko + Qo) + jf^oQoQo + jf^2{QiQi + Q2Q2). (4.5) 

As in the previous section, the Hartree approximation of any generator, the corresponding 
classical variable, will be denoted by the same symbol without a hat, the map of a product 
will be the product of the maps, but the map of a commutator will be i times the Poisson 
bracket of the two maps in the same order. We can now state the theorem that generalizes 
the one for n = 1. 

Theorem: For a class of many-Fermion Hamiltonians belonging to the enveloping algebra of 
Sp(2,R), of which Eq. (|4.5| ) is a prototype, there exists a six-dimensional family of Slater 
determinants, defined by density matrices, p(x, x'lQo, Pq, Qi, Pi, Q2, P2), that describe states 
belonging to an irreducible representation, which 

(i) induce a mapping of Sp(2,R) onto a symplectic manifold {Qq, Pq, Qi, Pi, Q2, P2) = {Q, P) 
in which the generators, G, are mapped as classical dynamical variables, i. e. , 

G^G{Q,P). (4.6) 

The mapping of products and commutators is as previously specified. 

(ii) The density matrices have the diagonal form 

p{x, x'\Q, P) = Y. MAQ, PWhiAQ, P), (4.7) 

h 

where the sum is over occupied orbitals, h. Furthermore il)h{x\QP) has the form 

MAQ^ P) = exp[tSix, P)]Mx, Q), (4.8) 
S{x, P) = \{xl + xI)Pq + \{xl- xl)Pi 

+X1X2P2, (4.9) 
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and (j)h is the solution of a constrained Hartree equation that will be specified and solved 
in the course of the demonstration. Altogether, these results will establish that the density 
matrix ( [4.7| ) is a solution of the time- dependent Hartree equation in the density-matrix form. 

We turn to the proof of this theorem. Since the technique is a relatively straightforward 
generalization of the proof given in the previous section, we shall be sparing of details, except 
in so far as these bring in something novel compared to the previous case. The first step is to 
replace the commutator algebra by a Poisson-bracket (PB) algebra. One new aspect is that 
we have an angular- momentum generator, L12. The PB relations between the coordinates 
and this quantity determine it to have the (not surprising) value 

Lu = QiP2-Q2Pi. (4.10) 

Next we consider the PB relations of the coordinates with the "deformation" generators, 
Sij. We thereby obtain nine differential equations which have the solutions, 

Su = {Pi + Po){Qi + Qo) + P2Q2, (4.11) 
^22 = (Pi - Po)iQi - Qo) + P2Q2, (4.12) 

512 = ^0^ + ^2^0- (4.13) 

A simple check on these results is that the sum of 5*11 and S22 must be a rotational invariant. 
We find 

Su + S22 = 2P0Q0 + 2PiQi + 2P2Q2. (4.14) 

With the above results, we can determine the momentum generators, Kij, by considering 
the nine PB relations of the coordinates with these generators, which are determined by the 
known values for Sij and L12. We thus find 

Pu = (Qo + Qi){Po + Pif + (Qo - Qi)P2 

+2Q2P2iPo + Pi)+XniQ), (4.15) 

P22 = (Qo - Qi)iPo - Pif + (Qo + Qi)P| 

+2Q2P2(Po - Pi) + X22(Q), (4.16) 

P12 = 2Q0P0P2 + 2Q1P1P2 + Q2(Po - Pi + Pi) 

+Xi2(Q). (4.17) 

Here, XijiQ) ^^^6 three unknown functions of Q that remain to be determined for a completion 
of the mapping. Notice in passing that the classical kinetic energy, 

K = i(Pn + P22) = QoiPo + Pi + Pi) 

+2Po(QiPi + Q2P2) + \{xii + X22), (4.18) 

is rotationally invariant, provided the sum of the last two terms on the right hand side has 
this property. 

For the determination of the unknown functions, Xij) we study the PB relations of the 
deformation generators Sij with the momentum generators Kij. As an example, the three 
relations involving Sij and P12 yield the differential equations. 
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(ni I ^ \i^X\2 . dxi2. . ^ dxi2 //I 1n^ 

+ ^'^^^ +907^ + = ^'-''^ 

^ dXi2 ^ dxi2 I . N 
where x = ^Xii + X22)- These equations have the solutions 



^ Q2 / \ 



whereas the remaining differential equations are compatible with the above and in addition 
yield 



X- = Uxn - X22) = g2 -Ql- Ql ■ ^^-2^) 



In order, finally, to determine the constant C, we calculate the map, C2, of the second 
order Casimir invariant in two ways. On the one hand we substitute the maps of the 
individual generators into the appropriate expression and thus find 

C2 = 2C. (4.25) 

On the other hand, we calculate the Hartree value of the invariant directly, as done in 
Appendix A. We thereby obtain (for a special representation) 

C = '^N\ (4.26) 

We complete part (i) of this demonstration by displaying the collective Hamiltonian that 
emerges from these considerations, 

He = QoiPo + Pi + P2) + 2Po{QiPi + Q2P2) 

+V{Q), (4.27) 



V{Q) = go + l/^oQo + \^2{Ql + Ql) + 



Qo~ Q\~ Q2 

(4.28) 



For part (ii) of the theorem, we must construct the density matrix that solves the time- 
dependent Hartree equation. The first step is to separate iIj{x\Q,P) into two factors, as 
done in Eq. ([4.8|) . From the time derivatives of the coordinates, calculated in two equivalent 
ways, from the classical equations of motion and from the Hartree averages of the quantum 
equations of motion, we obtain three equations, one for each coordinate, that generalize the 
single equation obtained by combining ( p.21| ) and (|3.22 ). These yield the solution 

S{x, P) = \{xl + xl)Po + Uxl- xl)Pi + X1X2P2. (4.29) 
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In analogy with Eqs. (|3.24| ) and ( |3.25| ), we introduce the orbital that is the solution to 
the time-dependent Hartree equation, iph, leading to the equation 

e,^, + ^ E (^M?^ + ^.^)=^V^'., (4-30) 

where the time derivatives are to be replaced by the classical equations of motion and the 
Hartree Hamiltonian, 7i, has the form 

H = ko + qo + KoqoQo + /«2(gi<5i + q2Q2)- (4.31) 

As before terms of second order in the classical momentum variables cancel. The zero-order 
terms yield the constrained Hartree equation, 

eh(ph = - E >^^q^)<Ph, (4.32) 
dV 

The first-order terms yield three equations that state, in analogy with Eq. ( |3.28| ), the iden- 
tity, when acting on between certain linear operators in the collective coordinates and 
corresponding linear operators in the single-particle coordinates. The statement and proof 
of these identities is given in Appendix C in order not to interfere with the flow of the main 
argument. 

We turn then to the solution of the constrained Hartree equation, ( [4.32| ). The operator 
that appears on the right hand side of this equation will be called 7Y. The cross-terms in the 
potential energy can be eliminated by an orthogonal transformation to the intrinsic system, 

(4.34) 
(4.35) 
(4.36) 
(4.37) 

(4.38) 

(4.39) 

(4.40) 
(4.41) 



Xl = 


cos^xi + sin^X2, 


X2 = 


— sintei + cos6'x2. 


COS 2^ = 


-{Col-ul)/Du 


sin 26 = 


-4A12/A, 


ul = 


C CQoiQo-Qi) 
Do Dl 


ujI = 


C CQo{Qo + Qi) 
Do Dl 


^12 = 


CQ0Q2 


Do = 


Qo ^ Qi ~ Q21 


Di = 


\JlQ\l, + {ul-ulY\ 



(4.42) 

These equations transform Ti to the form 

H = \{pl+ ujIxI) + \{pl+ ujIxI), (4.43) 

where 
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cUi2 = |C|/(Qo + gf, (4.44) 
^l = \C\/{Q,-Q)\ (4.45) 



Q = ^Ql + Ql = Qi- (4.46) 

The last of these expressions shows that the intrinsic system is the one in which Q2 vanishes. 

In consequence of the above transformation, the solution of the constrained Hartree 
equation is 

0,(x,g) = (cuic.2)^/VrHv^^i)0i?°Hv^^2). (4.47) 
Together with the material relegated to Appendix C this completes the proof of the theorem. 



B. Application of the theory of large amplitude collective motion 

Following the procedure outlined for the one-dimensional case, the first step is to compute 
the collective potential energy. Because in the intrinsic frame the Hartree Hamiltonian is 
a sum of two harmonic-oscillator contributions (with different frequencies), the procedure, 
including the use of the virial theorem, for obtaining a differential equation for the potential 
energy follows through without a hitch. Here it yields the partial differential equation 

V = 2Qo + -Ml + K,Q\) - - Qi^, (4.48) 



with the solution 



V = Q, + \Ml + K,Ql) + ^1^%. (4.49) 



In a general coordinate system, we should replace Qi by y <5i + Qi- Thus we recognize also 
in the present context that ^ is a scalar. 

Consider next the kinetic energy, in the form 

T = ^B^-P^P,, (4.50) 

where the indices take the values 0, 1, 2. We illustrate the two steps that enter into the cal- 
culation of the mass matrix elements B^^ . As one of three equations for the time derivatives 
of the coordinates, calculated in two ways, we have 

Qi = Su -822 = ^ = B^^P^. (4.51) 

As one of six equations necessary to complete the calculation, we have 

[Qi, 5n - S2.] = = 2Q, = B^\ (4.52) 

that combines a known PB relation with a derivative of ( |4.51| ). The completion of the 
procedure just exemplified yields the mass matrix previously determined by the theorem, 
that we display in matrix form. 
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B 



/2Qo 2Qi 2Q2 

2Qi 2Qo 
V 2Q2 2Qo 



(4.53) 



As explained in the previous section, this estabhshes once more that we have an exactly 
decoupled manifold. Again the alternative calculation within the framework of our theory 
of collective motion is summarized in Appendix B. 



V. MONOPOLE PLUS QUADRUPOLE MODEL ASSOCIATED WITH THE 

ALGEBRA OF SP(3,R) 

A. Classical mapping and Hartree solution 

In this section, that indeed quotes the results of most future interest to us, we shall drop 
the pretensions of formality adhered to in the previous sections. Except for very few points 
emphasized below, the results to be proved as well as the techniques used to carry out the 
demonstrations should be evident by now. We divide the generators into monopole and 
quadrupole parts. In order to keep better track of the significance of the variables, we adopt 
an alphanumeric subscript notation, illustrated by means of the one-particle operators. 



13 = :^{xl + xl + xl), (5.1) 

qdi^-l{xl + xl-2xl), (5.2) 

(ld2 = ^{xl-xl), (5.3) 

qoi = X1X2, (5.4) 

qo2 = X1X3, (5.5) 

Qos = X2X3, (5.6) 

Ss = :^{xiPi + X2P2 + X3P3), (5.7) 

Sdl = ~{xiPi + X2P2 - 2X3P3), (5.8) 
etc., 

ks^lipl+pl+pl), (5.9) 

kdi^-l{pl+pl-2pl), (5.10) 
etc. 



Associated with these one-particle operators are second quantized operators, designated 
by the same symbols in the upper case and carrying hats, and the Hartree maps of these, 
designated by the uppercase but hatless. By means of the PB algebra these are expressed 
as functions of six canonical pairs, {Qg, Ps), ■ • • , iQo3, -Pos)- The resulting espressions, which 
were calculated with the aid of the program, Mathematica, are too long to quote in their 
raw form. To obtain more concise expressions, we revert to a spherical tensor notation 
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Qo = 


Vs. 


(5.11) 


4 = 


7KiQo2+iQo3), 


(5.12) 




—7^{Qo2 — iQos), 


(5.13) 


Q2 = 


-^{Qd2 + iQoi), 


(5.14) 




-^{Qd2 - iQoi), 


(5.15) 



for Q, and similar definitions for S and K. In the following expressions, the classical limit 
of Qm is denoted by S^, and the conjugate momenta are denoted by 11^. We use Qs and Ps 
for the monopole coordinate and its momentum. The spherical components of -S" are denoted 
by E, etc.. Finally, Lm will refer below to the standard spherical components of the angular 
momentum vector. 



We thus find that 



Xs 
C2 



iVlO [E X n] 



(1) 

m ' 



-{E-U + QsPs), 



4 



9 

V5 



15 
6\/5 



[n X Uf^ X E 



(2) 



[n X n] 



(2) 



X 



(2) 
m 



[nxn]^^) xs' 



(2) 



= -Cdg^ In (4Q3 - 4S • Sg, - 73(2)) , 
= In {AQl -AE-EQ,- I,{E)) , 

= -30sgs-2^$^5^ = 9C, 



(5.16) 

(5.17) 

(5.18) 



(5.19) 



(5.20) 
(5.21) 
(5.22) 

(5.23) 
(5.24) 



He = -{K, + Qs) + \kqQI + \k2E ■ E + \klL ■ L. 



(5.25) 



The above expressions contain the standard notation for angular momentum coupling. We 
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recognize Xs and Xm as the singular functions that arise in the determination of the momen- 
tum generators. Here the partial signs represent derivatives with respect to the variables 
written, and C is the integration constant that has appeared analogously in the calculations 
for lower dimensionality. A complication that has not appeared previously is the occurrence 
of the cubic invariant, J3, constructed from the tensor S^- Finally, from the given form 
of the collective Hamiltonian, He-, it is obvious what starting many-body Hamiltonian was 
used. Compared to the corresponding two-dimensional case, we have only added an angular 
momentum coupling. 

To construct the density matrix, we write 



i)h{x\Q, P) = exp[iS{x, P)](t>h{x, Q), 



and find 



s{x, p) = PsQs + n ■ g. 

For the real orbitals, we find 



i=l 



(5.26) 



(5.27) 



(5.28) 



as the solution of a constrained Hartree equation, referred to the intrinsic or barred system 
(Cf. (|4.44| )). For the constrained Hamiltonian we write 



dV 

Ws 



dV 
Wo 



+ qd2 (/«2<5 



dV 



d2 



OQ 



d2 



(5.29) 
(5.30) 



1 + -^{i^gQs 



dV 



1 



dQs 
dV 



2 

1 + 7;{i^oQs 



--^('«2<5o 
2 

2 

+ -^{l^2Qo 



dQo 
dV 



dQs 
dV 



dQo 
dV 



dQs 
dV 



dQ, 



+ (K2Q, 



dV 



d2 



dQ 



- {l^2Q 



d2 



dV 



d2 



dQ 



d2 



(5.31) 



(5.32) 



(5.33) 



Alternative expressions for the frequencies that may be found by substituting the derivatives 
of the potential energy will not be recorded here. 

In addition to the above relations and conditions, there are also equations that relate 
certain linear differential operators with respect to the collective coordinates to linear dif- 
ferential operators with respect to the intrinsic coordinates. These were discussed for the 
two-dimensional case in Appendix C, but will not be discussed at all for this case. 
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B. Application of the theory of large amplitude collective motion 



Following the procedure developed for the one- dimensional case and previously also ap- 
plied to the two-dimensional case, we obtain a partial differential equation for the potential 
energy. In stating this equation and its solution we shall suppress the bar notation, under- 
standing that the collective coordinates refer to the intrinsic system. We thus find 

This equation has the solution 

V=IQs + l^oQl + |'^2(g^ + Ql) + Xs, (5.35) 



where Xs is the singular scalar function displayed in Eq. ( |5.21 ). In the present calculation 



the quadratic and cubic scalars, I2 and I3 are found in the versions to which they reduce in 
the intrinsic system, namely, 

l2 = Ql + Ql, (5-36) 
iQl-^QoQl)- (5.37) 



3v^^ 

There remains only the problem of computing the mass tensor. This remains a simple 
algebraic task, if one applies the method first described for the one-dimensional case. The 
calculation is simplest to carry out in terms of Cartesian variables, but the equivalence of the 
result to that given by Eqs. (5.19) and (5.25) can then be verified. Nothing new is learned 
by repeating either the calculation or the results. 



VI. CONCLUDING REMARKS 

In this paper we have developed a new exactly solvable example for the theory of large 
amplitude collective motion by looking at the classical limit of the algebra of Sp(3,R) and 
a Hamiltonian defined within the enveloping algebra of the algebra. Possibilities exist both 
for the application and for theoretical refinement of the results of this paper. A natural 
first application would involve the requantization of the classical collective Hamiltonian and 
comparison of the results of diagonalizing the resultant Bohr Hamiltonian with the results of 
an exact diagonalization of the corresponding many-particle Hamiltonian |]l4l . A theoretical 



refinement would be to attempt to upgrade the classical mapping to a full quantum boson 



mapping |T^. This is probably not difficult for the one- and two-dimensional cases, but 
might prove laborious for the interesting n = 3 case. 

It would be inappropriate to conclude this paper without mentioning possible connections 
with other work. For instance there has been extensive research on the unitary representa- 
tions of the non-compact symplectic algebras considered in this paper [[IGHl^, that can be 



viewed as exact boson mappings of these algebras. The relation of our special mapping to 
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this work may be worth pursuing. Along a different hne, we may be said to have produced 
sohton solutions to a highly simplified class of field theories, using Hamiltonian methods. 
It would be interesting to investigate if our methods can teach us something about the 
Hamiltonian approach to other soliton models ||19| . 



This work was supported in part by the U. S. Department of Energy under Grant No. 
40264-5-25351. 



APPENDIX A: EVALUATION OF SECOND ORDER CASIMIR INVARIANT 

Consider the one-dimensional case. The evaluation is based on the idea that since the 
monopole-monopole interaction belongs to the enveloping algebra, changing its strength 
does not change the irreducible representation. In particular, we may set the strength to 
zero. In that case the ground state, that defines the representation of interest to us, is the 
Slater determinant obtained by filling the lowest orbitals. Since in this state the Hartree 
value of S vanishes, we have from Eq. (2.29) 



N-l 

C2(irf) = -(E^/^)' = -(E 

h 

^ -N\ (Al) 
4 

In this evaluation, as in the ones that occur below, it is strictly the definition of the clas- 
sical variables as Hartree averages that intervene. Furthermore, we have utilized the virial 
theorem for the harmonic oscillator described by the Hamiltonian ^{p"^ + x"^) . 

For the two and three-dimensional cases, we shall evaluate the Casimir invariant only 
for the closed-shell and for the representation containing the ground state. This allows us 
to take fullest advantage of circular and spherical symmetry, respectively. We provide a few 
details for three-dimensions only. We have 

3 
i=l 



For particles, this becomes 



1 



(6Af)i/3 



C2(3rf) = --( E I 







= _1qV^NS/3_ (A3) 
16 ^ ^ 

Here we have used the fact that the level with energy n has degeneracy ^n"^ and that the 
number of shells needed for particles is (6A^)^/^. These, of course, are approximate values 
needed to get the answer to leading order in A^. 

The corresponding calculation for two dimensions yields 
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C,(2d) = -^N\ (A4) 

In future applications, we shall be interested in irreducible representations associated with 
open-shell nuclei. The generalization to such cases of the elementary calculations just pre- 
sented will be deferred to the occassion when those applications are presented. 



APPENDIX B: ALTERNATIVE COMPUTATION OF THE MASS TENSOR 

In the body of the text we have described a new method for computing the mass tensor of 
a decoupled or approximately decoupled manifold. This method is particularly convenient 
when the cranking operator or operators can be represented in a basis of single-particle 
operators that depend only on the coordinates, though it can also be applied with some- 
what increased effort in more general cases. Nevertheless, it is important to show how the 
same results can be obtained by the methods developed in our previous work. In fact, the 
monopole model was instrumental in leading to a necessary generalization of the theory of 
the mass tensor as it had been presented and utilized by us. The generalization required 
is presented in some detail in 0. We do not wish to repeat the technical details here, but 
rather we shall "remind" the reader of the basic idea and be content with the quotation and 
application of the final result. 

The theory of large amplitude collective motion is based on the idea that starting from 
a Hamiltonian, quadratic in the momenta but otherwise arbitrary, one can introduce a 
point transformation chosen so that in the new coordinates the existence of a decoupled 
coordinate manifold is either apparent, or in more realistic cases can be demonstrated to be 
approximately true. The restriction to a coordinate or point transformation was based on the 
observation that this type of transformation maintained exactly the quadratic dependence 
on the momenta. However, the formulas derived from this procedure failed to reproduce the 
correct mass in the monopole example. Thereafter, it was discovered |^ that for consistency 
to second order in the momentum one had to replace point transformations by more general 
canonical transformations, correct to second order in the momenta. In the paper cited 
details were supplied only for the case of one collective coordinate, but for the separable 
Hamiltonians considered in the present paper it is straightforward to generalize the result 
to any number of collective coordinates. Below we describe the result of this generalization. 

We suppose that there is a set of self-consistent cranking operators, 

Q^ = tr(gV). (Bl) 
It follows then that the component B^^ of the mass tensor is given by the formula 

Remarkably, where the present formula contains the cranked Hartree Hamiltonian, 7i, the 
previous, incorrect formula contains just the unconstrained Hartree Hamiltonian Ti. Pre- 
cisely because the single-particle orbits referred to in the above formula are the eigenmodes 
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of the constrained operator, and using completeness as well, the formula for the mass tensor 
can be transformed into the simple formula involving double commutators. 



B'^ = \ii{p[[q\nlq^]+i^3}. (B4) 

For the cases treated in the text, the are coordinate operators. Consequently 71 may 
be replaced by the single-particle kinetic energy. It is then a trivial calculation to show that 
all the results of the text are reproduced. 



APPENDIX C: SOME ADDITIONAL DETAILS CONCERNING THE SOLUTION 
OF THE TIME-DEPENDENT HARTREE EQUATION 

In the construction of the density matrix that solves the time-dependent Hartree equation 
for the two-dimensional case, we bypassed the study of a set of conditions that must be 
satisfied by the real orbitals (ph in addition to the constrained Hartree equation ( [4.321) . 



These conditions, here three in number, arise from terms linear in the collective momenta 
in the time-dependent Hartree equation. They are of the form, 

on , on ^'^ ^nQ^ 



To show that each orbital, (ph, given by Eq. (f4.32|) satisfies these conditions, it is conve- 



nient to transform them to the intrinsic system. In addition to the transformation equations 
( [4.34| ), ■ ■ -, we introduce polar coordinates 

Qi = gcose, (C4) 
Q2 = Qsme, (C5) 

that are recognized as intrinsic coordinates, if we make the identifications Q = Qi and 
9 = —29. Straightforward calculation then permits us to replace Eqs. (PT|)- ( |C3D by the 
equations 

on , on - dcp d4> 

2Qo 50 _ - 50 90 

-Q-^e-''''W^~'''W, ^^^^ 



To verify the last set of equations, we utilize the explicit form ( 4.48 ) of the orbitals. 



as well as the Eqs. (|4.45|) : and (|4.46| ) that relate the barred frequencies to the collective 



coordinates. We must also remember that the barred coordinates are also functions of the 
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collective coordinates, according to the equations that are the inverses of ( [4.34| ) and (f4.35|) 



and the relation between 9 and B. After some algebra, we find that we can duplicate (|C6| 
and ( |C7| ), but in place of (|C^), we deduce the result 



2— = -xi— + X2 — , (C9) 
dQ dx2 dxi ' 

which obviously describes correctly the rotational properties of an orbital, in view of the 
connection between and 6. To verify ( |C8|) , finally, we calculate the "commutator" of Eqs. 
and ( p9D and indeed deduce the desired result. 
The corresponding calculations for three dimensions will not be reproduced here. 
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